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SUMMARY 


An analytical investigation has been made to establish the 
influence that different mathematical and aerodynamic models have 
on the computed spin motion and to determine the importance of 
some of the aerodynamic and non-aerodynamic quantities defined in 
these models. Because of the knowledge acquired in this investi- 
gation, it should be possible to identify the reason for any dis- 
crepancy that may be noted in the future between predicted and 
flight recorded spin motions. 

« 

The conventional analytical technique used in spin investi- 
gations does not include the aerodynamic forces and moments acting 
on a spinning aircraft due to steady rotational flow nor does it 
limit the contribution of the rotary derivatives to the oscillatory 
component of the total angular rates. Whereas, an analytical tech- 
nique employing rotation-balance data investigated herein does both. 

It was shown that a spin cannot be computed using the conventional 

^ical technique when a rotary derivative is unstable over some 
portion of the angle-of-attack region. This was not the case when 
the analytical technique based on rotation-balance data was used. 

Some of the conclusions arrived at in the past using the conventional 
analytical technique, therefore, are incorrect relative to the 
significance of rotary derivatives on the spinning motion. Also 
the spin computed with rotation-balance data duplicated the developed 
spin obtained in the spin tunnel whereas, an appreciably less severe 
spin was realized with the conventional technique. It was indicated, 
therefore, that the aerodynamic moments generated in a spin due to 
rotational flow, as measured by a rotation-balance, are indeed signi- 
ficant. The analytical technique used to date is, therefore, a 
pseudo-technique which cannot be used to predict aircraft spins. 

A preliminary study indicated that static aerodynamic derivatives 
can be extracted from rotation-balance data. It may be possible, 
therefore, to simply add the aerodynamics associated with steady 
rotation flow to the conventional aerodynamic model. 

It was also shown that during experimental-analytical correla- 
tion studies the flight recorded control time histories must be 
faithfully duplicated since the spinning motion can be acutely sen- 
sitive to a small change in the application of the spin entry con- 
trols. However, an error in the assumed inertias, yawing moments at 
high angle of attack and initial spin entry bank angle do not influ- 
ence the developed spin significantly. The damping in pitch deri- 
vative and center of gravity location were shown to play a signifi- 
cant role in the spinning motion. ° 

It was also concluded that the experimental spin investigations 
conducted in a constant atmospheric density environment duplicate 
th§ Froude number only at the initial full-scale spin altitude since 
the full-scale airplane at high altitudes experiences large density 
changes during the spin. Therefore, the full-scale rate of sink and 
spin rate predicted from a constant density environment model test 
will be greater (conservative) than the actual airplane values. 


* 
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INTRODUCTION 


The Langley Research Center has a broad research program 
designed to advance the state-of-the art in the area of stall/ 
spin technology. One major requirement existing in this area 
is the development and validation of reliable theoretical 
methods for prediction and analysis of stall/spin characteris- 
tics. Although theoretical studies have been made in the past, 
no concentrated nor continuing effort has been made to rigorously 
validate the analytical techniques employed in these studies. 

In view of the urgent need for valid theoretical techniques, this 
study was performed as a step toward advances in this area. 

The analytical technique most commonly employed relative to 
the study of the stall/spin phenomenon Involves the simultaneous 
solution of the equations of motion and associated formulas. How- 
ever, little confidence is placed In the technique and this atti- 
tude will persist until the incipient, developed and recovery 
phases of the spinning motion recorded in flight with the full- 
scale airplane are truly matched by a computer solution. 

. The inability to satisfactorily demonstrate this "motion 
matching" in the past might be attributed to the use of an incom- 
plete or a poor representation of the stalled aerodynamics that 
exist in a spin. For instance, the aerodynamic model employed 
analytically may be based on low Reynolds number wind tunnel tests 
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when it is known that Reynolds number has a significant effect on 
the aerodynamic characteristics obtained at spinning conditions 


for some current configurations. Also, the technique employed 
to represent the aerodynamic forces and moments may very well 
be a pseudo-technique which must be abandoned when attempting 
to compute spins in which the aircraft rapidly rotates about a 
vertical axis located near the center of gravity. For these 
spins, the sideslip angle varies along the length of the air- 
craft; being of opposite sign forward and aft of the axis of 
rotation. It may be necessary, therefore, that a rotation- 
balance rig be employed to obtain a set of aerodynamic data 
while the model is under the same local flow conditions that 
exist during the actual spin. 

The investigation reported herein is the first phase of a 
two part effort. The objective of the final phase is to deter- 
mine how well analytically determined motions correlate with 
those measured in flight when the theoretically and experimen— 
tal ly determined motions are based on the same Reynolds number. 

This is to be accomplished by using i — scale radio— control 
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model time histories. To be able to quickly identify the cause 
for any discrepancy that may be noted during this forthcoming 
correlation study, it is necessary that the importance various 
factors play in the computed motion be known. Therefore, the 
objective of the first phase of the investigation is to establish 
the required background of knowledge. 

TECHNICAL APPROACH 
Required Studies 


As indicated in the Introduction, the principal objective of 


this investigation is to establish a background of knowledge 
relative to 

1) computational procedures and analytical models (both 
mathematical and aerodynamic) and 

2) the influence that various non— aerodynamic and aero- 
dynamic quantities have on the computed spin motion 

to aid in identifying the reasons for a discrepancy between 
predicted and flight spin motions. Also, by achieving the ob- 
jective one can select an analytical model which appears to give 
solutions more representative of the "real world". The following 
four distinct studies were required to achieve the objective. 

The first study attempts to verify that the computational 
procedures employed would not in themselves be the source of a 
discrepancy between computed and measured motions. The concern 
is that an undetected error may have inadvertently been incorp- 
orated into the large, complex computer program required for cal- 
culating spin motions or that the integration routine employed 
would be lacking relative to the rapid motions experienced with a 
scaled model. 

The second study determines the significance of not duplica- 
ting the attitude of the aircraft in space at the beginning of the 
spin entry maneuver or the control time histories. The first item 
is of concern since attitude gyros are not included in the instru- 
mentation package placed aboard Langley radio— control models . The 
latter item is to be studied since instrumentation errors, aero- 
elastic effects, engineering approximations of complex recorded 
flight control time histories, etc., may preclude simulating the 
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exact control time histories experienced by the aircraft. 

The third study determines the influence that possible errors 
in the physical (i.e., moments of inertia, center— of— gravity lo- 
cation) and aerodynamic characteristics selected to represent the 
aircraft have on the computed motion. 

These three studies used the conventional analytical technique 
where the aerodynamics are represented by static data and dynamic 
derivatives. The aerodynamic model used for these studies is re- 
ferred to as model 1 which is described under the appropriate head- 
ing. The study concerned with the influence of aerodynamic charac- 
teristics also investigated a set of aerodynamic data referred to as 
model 2. 

The fourth study conducted herein investigates the signi- 
ficance of employing a different type of mathematical and aerody- 
namic model which is described in reference 1. In this instance, 
the models are based on available rotation-balance data. 

This technique is investigated since the conventional technique 
does not include the aerodynamic forces and moments acting on the 
spinning aircraft due to steady rotational flow nor does it res- 
trict the dynamic derivatives to the oscillatory component of the 
total angular rates. Whereas, the technique of this latter study 
does both (see Appendix A), The aerodynamic models employed dur- 
ing this study are referred to as models 3 and 4 which are also 
described in this section. 

Analytical Technique 

The influence that a mathematical or aerodynamic model, or 
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some individual non-aerodynamic or aerodynamic quantity has on the 
computed motion was determined by examining the path of the center 
of gravity and the motion about the center of gravity during the 
incipient and developed spin phases. All of the studies referred 
to herein, therefore, required the computation of large angle mo- 
tions using a large angle, six-degree-of-freedom, plotted output 
digital computer program with nonlinear, multifunctioned aerodynamics 
coefficients. This program solved the equations of motion and 
associated formulas presented in Appendix A. 

Forty second time histories, deemed of interest to the reader, 
are presented herein. One 120 second time history is also included 
for calculation made to simulate the ful-scale aircraft. Each time 
history is presented in a figure consisting of two 11” x 17" pages 
and each page contains the following variables plotted versus time. 

Page 1 Page 2 

r p 

cc q 

y3 r 

V R Si. 

It was felt that the spinning motion was sufficiently illus- 
trated by describing the relative wind vector at the center of 
gravity on the first page and the angular velocities about the 
axes having their origin at the center of gravity on the second 
page. 

Aerodynamic, Inertia and Mass Data 

The aerodynamic data which represented the clean, 22 degree 
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leading edge wing-sweep configuration of a - - scale 
radio-control model are presented in Appendix C. All the data 
in this appendix were reduced to standard coefficient form on 
ohe basis of the following geometric characteristics: 

s = 5.65 ft 2 c = II.76 in b = 6.4l ft 

and were transferred during the computations to a 14 percent c 
moment reference center to correspond with the flight model 
center-of-gravity location. 

The Appendix C data were employed to construct four basic 
aerodynamic models which, as previously mentioned, represent two 
different type of models i.e., the conventional aerodynamic model 
and one based on rotation-balance data. 

C The conventional models are identified in Appendix C as models 

1 

1 and 2. Model 1 has been used in a previous investigation (refer- 
ence 2) and includes static data, experimentally obtained over a 
wide range of Reynolds numbers in various test facilities, as well 
as estimated rotary derivatives. These data are referred to a 
moment reference center located longitudinally at 9 percent C unless 
otherwise noted. Model 2 consists of static data (some of which is 
presented in reference 3) and rotary derivatives (reference 4) 
which were experimentally obtained in the Langley Full-Scale Tunnel 
facility with a - scale model. As shown, rotary derivatives 
were measured at two forced frequencies. However, for this investi- 
gation, the rotary derivatives were assumed to be invariant with 
frequency and the values obtained at the high frequency were selected 
as the base for model 2. The model 2 data are referred to a moment 
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center located longitudinally at 16.25 and 31.55 percent c for 
the static and dynamic data, respectively. All of these data 
were measured at a Reynolds number based on c between . 4 and .5 x 10 ^ 
which is equivalent to the flight Reynolds number of the radio- 
control model. 

Models 3 and 4 both use the same rotation-balance data but 
different rotary derivatives. Model 3 employs the estimated 
rotary derivatives from model 1 whereas model 4 used the experi- 
mentally derived derivatives from model 2 . The rotation-balance 
data were obtained at a Reynolds number based on c of .47 x 10^ in 
both the clockwise and counterclockwise direction on the rig shown 
in figure 1 for the following control configurations. 

all neutral controls (i_, t> . S = 0°) 

s a R ' 

stick full aft (i„ = -30°; & , = 0°) 

23 a K 

right pro-spin controls (i - - 30 °, 6 =4 5 °» 6 = - 30 °) 

s a R 

left pro-spin controls (i s = -30 ,-6 = - 5°* 6 = 30°) 

& R 

Data were recorded at 6 = 55 * 60, 65, 70 , 75 , 80, 85 and 90 
degrees for 0 values of 0 and 5 degrees at ..■ftk values of 0, . 0493 , 
.1030, .1523, .2015, .2553 and .3045 and were measured about a moment 
reference center located longitudinally at 31.55 c. These data 
are presented in Appendix C and were inputed "as is" into the com- 
puter in tabular form i.e., tables of aerodynamic data as a function 
of 0 and fl exist for each control configuration at 0 = 0 and 5 
degrees. The test attitude angles © and 0 were assumed to be equal 
to oc and , respectively when these tables were used. However, 
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since a = tan *” 1 (tan 0 e cos 0 e ) and y 3 = tan " 1 (sin © e sin 0 e ), 

this assumption results in some small errors when the 0=5 degree 
tables were used as shown below: 


0=5 


6 

r ✓ 

CL 

fi 

55 

54.90 

4.09 

60 

59.91 

4.33 

65 

64.92 

4.53 

70 

69.93 

4.70 

75 

74.95 

4.83 

80 

79.96 

4.92 

85 

84.95 

4.98 

90 

90.00 

5.00 


angle of attack and/or sideslip angle were computed 
momentarily which was greater or less then those available in the 
tables (overflow condition) the last value used would be maintained 
until the computed angle fell within the confines of the available 
data. 

A weight of 149 pounds, a center-of-gravity location of 14 
percent c and inertias of 2.205, 8.074, 10. 096 and 0 slug-ft 2 for 
X x> ly* ! z and I xz , respectively were experimentally determined for 

the i6 ** SCale radi °- contro1 model. (The corresponding full-scale 
values for the weight and inertias are 46,000 pounds, 72,700, 

266,000 and 333,000 slug— ft ). These model values were used, there- 
fore, to compute the spinning motion for a model whose mass is 
distributed slightly more in the fuselage than in the wings as 
indicated by the values of the inertia parameters, i.e.. 
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2 2 p 

( X x ~ I y)/ mb > (!y - I z )/mb and (I z - I x )/mb 2 equal -309, -106 
and 415 x 10~\ respectively. 

Initial Spin Entry Conditions 

All the studies conducted with the classical aerodynamic 
models used the following entry conditions unless otherwise noted. 
Motions were initiated at an altitude and speed of 4000 ft and 
173 ft/sec, respectively in a trimmed lg wings level attitude. 

The trim values for models 1 and 2 were as follows: 



Model 1 

Model 


-1.9° 

-0.29' 


6.1° 

O 

CO 

• 

I*- 

T 

17.5# 

19 M 


An accelerated stall entry was then performed. The entry 
maneuver involved applying a stabilizer deflection of -30 degrees 
and a differential horizontal tail deflection of +12 degrees at a 
rate of 9^.8 deg/sec starting at time equal to zero and 0.5 seconds, 
respectively. These maximum control deflections were then held con- 
stant during the ensuing motions. The rudder was left undeflected. 

A slight variation in this spin entry maneuver was employed in the 
aerodynamic study. In this instance, upon reaching a -30 degree 
stabilizer deflection the stabilizer was immediately reduced to 

-23 degrees at the 9^.8 deg/sec rate; and the differential horizontal 

* 

tail movement began at a time equal to 0.316 seconds rather than 0.5 
seconds. (It should be noted that a control rate Of 9^»8 deg/sec and 
a time of 0.316 seconds is equivalent to 30 deg/sec and 1 second. 
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respectively for the full-scale aircraft.) The thrust was 
arbritrarily reduced to zero above ?0 decrees angle of attack. 

The studies conducted for aerodynamic models 3 and 4, 
which are based on rotation-balance data, employed initial con- 
ditions which would be realized during the incipient spin phase. 
This approach was necessary since rotation-balance data were not 
available below an angle of attack of 55 degrees or for a side- 
slip angle greater than 5 degrees. The initial values at time 
equal zero seconds employed for these studies when spinning to 
the right were as follows: 



Model 3 

Model 4 

T 

0 

0 

V R 

93.4 

92.4 

H 

3939 

3912 

r 

- 61.6 

-63.9 

(X 

69.6 

72.6 

fi 

0 

0 

p 

0.7314 

1.1064 

q 

- 0.00045 

0.0345 

r 

2.0452 

2.4011 

0 

0 

0 

V 

-44.5 

-64.0 

Sa 

+5 

+5 


-30 

-30 


-30 

-30 


These initial conditions for models 3 and 4 were based on the 
values realized at t = 6 seconds during a 60 degree banked spin 
entry maneuver that had been computed with models 1 and 2, res- 
pectively. In retrospect, it appears that assumed initial values 
would have sufficed for this study. The time histories associated 
with this study are terminated before 40 seconds have elapsed. 

The reason being that computations were terminated upon reaching 
sea level. The control deflections at t = 0 seconds were predica- 
ted on the pro-spin control configuration for which rotation- 
balance data were available. Because of adverse yaw due to the lat- 
eral control, crossed controls are considered to be pro-spin for 
this model. 


RESULTS AND DISCUSSION 

Validity of Computational and Model Testing Technique 

As stated previously, the first objective of this investi- 
gation was to eliminate the possibility that an improperly program- 
med computer would be the source of some discrepancy that might be 
noted between computed and measured spins. A check case was selec- 
ted, therefore, which would thoroughly exercise a large-angle six- 
degree-of-freedom computer program. This check case was then solved 
using different programs and computers that were available in-house 
9-t Grumman and Langley . The solution from these facilities were 
compared and the studies presented herein were not initiated until 
perfect agreement was noted. 

Although the final objective of this investigation is limited 
to determining the degree of correlation obtained between computed 
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and measured model motions, there is always the underlying inter- 
est in how well the model motions simulate the full-scale aircraft. 
The time was taken, therefore, to briefly investigate the importance 

on the spin of a factor that is not simulated during experimental 
Langley spin studies. 


Dimensionally scaled Langley radio-controlled drop models and 
models tested in the free-spinning tunnel facility are dynamically 
ballasted for some full-scale altitude (usually 35,000 to 40,000 
ft) so that the models and full-scale airplane have similar geo- 
metrical paths; even though the radio-controlled models are tested 


near an altitude of 3,000 ft and the free-spinning models at sea 


level conditions. This is accomplished by making the weight of 

^ M T.T /*T\-3 


the model equal to 


tia equal to 'Jy 1 I A (N) 


~TK 


TT 


-5 


Wa (N) - - 5 and the model moments of iner- 


(N in the formulas represents the 


scale i.e . > i , i. etc.) In this manner the Froude number of the 
10 20 

full-scale airplane is duplicated i.e., similarity between inertia 
and gravitional forces is maintained. Consequently, the attitude 
angles ($ e , 0 e ,(/ e ), and turns required for recovery of the 

airplane and model are the same, the linear velocities of the air- 


plane center of gravity (c.g.) equals V M (N)“ 1//2 and the 


M 


v 1/2 


airplane angular velocities about the c.g. equalsn M (N)' 

It should be noted, however, that the Froude number is only 


duplicated at the initial full-scale spin altitude since the 

experimental spin investigation are conducted in a constnat 

♦ 

density environment whereas the full-scale airplane at high 
altitude experiences large density changes during the spin. 

The role that the rate of change in atmospheric density with 
altitude plays in the spinning motion of the aircraft, therefore. 
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is of interest. • 

The effect of density can be determined from figure 2 which 
presents the time histories ensuing when the same control manipu- 
lations are applied to a full-scale aircraft at 35 ,000 ft and 

a - - scale model at 4,000 ft. It can be seen that the over- 

10 

all incipient and developed phases of the full-scale spins are 
predicted rather well by the model. If one applies the previously 
discussed scaling factors, it will be found that V R , CC-, /3, p, q, 
r, and _n_ vs time are faithfully duplicated up to approximately 40 
seconds (full-scale value). 

Beyond 40 seconds the velocity traces between the model and 
full-scale aircraft depart. Since an aircraft decends through the 
atmosphere at a constant value of dynamic pressure, the resultant 
velocity (rate of descent) continuously changes and is an inverse 
function of rate of change of density with altitude. Consequently f 
the altitude loss and ratio of altitude loss per turn obtained with 
the scaled models becomes correspondingly inflated with decreasing 
altitude relative to the full-scale value. 

The rate of change in atmospheric density experienced by the 
full-scale aircraft at high altitude influences the a. and /) values 
in the developed spin only slightly. However, the angular veloci- 
ties are significantly affected. Since a flat spin is computed, this 
effect is most obvious in the yaw rate which shows the full-scale 
aircraft reaching an equilibrium value sooner and of a lower magni- 
tude than that realized with the model. The importance of not dupli- 
cating the Froude number during the developed spin phase is illus- 
trated by taking the model values at t = 31.6 sec (model time) and 
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scaling them up to the full-scale values which are then compared 
with the values computed at 100 sec for the full-scale aircraft. 


At t = 100 sec 

a 

A 

V R 

P 

q 

r 

JX 

Scaled-model values 

83 

-0.7 

303 

.300 

-.024 

2.69 

2.70 

Full-scale values 

81 

-1.6 

208 

.320 

-.003 

2.16 

2.20 


As shown, the V R and it of the developed spin measured in a constant 
density environment will be greater than the values realized with 
the full-scale aircraft and, therefore, are conservative. 

Importance of Duplicating Control Time Histories 

The sensitivity of the spinning motion to simple variations 
in the spin entry control time history was investigated by systema- 
tically changing the time of control application (t Q ), the control 
rate or control deflection of the stabilizer and differential hori- 
zontal tail. It was found that although the incipient spin motion 
reflected these control changes, as would be expected, the developed 
spin remained unaffected until some small additional control change 
triggered a radically different type of spinning motion. This 
ex'treme control sensitivity is illustrated by the time histories 
presented in figure 3 and summarized for various control manipu- 
lations on the next page. 
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Fig. Control Manipulation Ensuing Motion 



*0 



t 

o 

B a 

6 

a 



sec 

deg 

deg/sec 

sec 

deg 

deg/sec 


3a 

0 

-23 

94.8 

1.0 

+12 

94.8 

post stall gyration 

3b 

0 

-23 

94.8 

.316 

+12 

94.8 

steady flat spin 

3c 

0 

-30 

30 

.316 

+12 

94.8 

moderate osc.spin 

3d 

0 

-30 

30 

.500 

+12 

94.8 

steady flat spin 

3e 

0 

-30 

94.8 

.500 

+ 5 

94.8 

steep osc.spin 

3f 

0 

-30 

94.8 

.750 

+ 5 

94.8 

steady flat spin 


It can be seen from the time histories and the above table 
that an apparently insignificant change in the control time history 
can alter radically the nature of the ensuing motion. The reader 
is reminded that this set of data was specifically selected to 
illustrate this point and that actually broad bands of control 
manipulation exist for which the developed spin remains unaltered. 

The boundaries separating these broad bands ( and correspondingly 
different motions) are extremely narrow. The boundaries being 
defined by some critical control phasing. Unfortunately, the con- 
trol phasing required to trigger a radical change in the behavior 
of the airplane cannot be predicted at this time. For example, 
the steep oscillatory spin obtained after the application of + 5° 
of differential horizontal tail was not suprising. The encounter 
of a flat spin by delaying the application of this small control 
deflection by 0.25 seconds obviously was unanticipated. 

In light of these results, it is recommended that analytical- 
experimental correlation studies employ flight time histories which 
•are specifically obtained for such studies and that considerable 


thought be given to the type of control manipulation used for spin 
entry. Also, if the nature of the analytical spin is appreciably 
different from that obtained in flight, the effect of control per- 
turbations about the flight measured control manipulation should be 
investigated analytically. 

Importance of Duplicating Initial Bank Angles 

The influence on the spinning motion of having the aircraft 
banked at the time spin entry controls are applied was determined 
for initial bank angles of 0, + 10, + 40, and + 60 degrees. It was 
found that the developed spin characteristics are unaffected by the 
initial bank angle, but that a noticeable influence on the incipient 
spin motion is realized for bank angles 40 degrees and greater. The 
effect on the incipient spin motion can be seen by comparing the time 
histories presented in figures 4a and b (0 = 6 0 ° and - 60 ° respec- 
tively) with figure 3d (0 O = 0°). The 0 Q = -60° does not differ 
significantly from the 0 Q = 0° spin entry whereas the 0 O - 60 ° exhi- 
bits a more rapid decay in the c*, /3 , p, q, and r oscillations. 

Entry into the developed steady flat spin, therefore, is achieved more 
smoothly by having the aircraft initially banked in the direction 
of the spin. It is interesting to note that this is exactly the 
situation that is experienced by a pilot when he picks up a wing 
with a lateral control that generates adverse yaw. 

Significance of Errors in Assumed Physical Characteristics 

Sensitivity of the spinning motion to center-of -gravity 
location can be determined by comparing the time histories presented 
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In figures 5a, b and 3d for a c.g. location of 6, 26 and 14 
percent c, respectively. It is shown that an aircraft capable 
of entering a flat spin will do so regardless of the magnitude 
of the spring constant. Of course, as would be expected, there 
is a small decrease and increase in the spin rate of rotation 
and angle of attack, respectively as the c.g. moves aft. The 
most significant effect, however, is on the incipient spin phase of 
the motion. As shown in figure 5, the effect of c.g. location on 
the time required and the type of motion experienced before the 
equilibrium spin condition is reached is appreciable. Moving the 
c.g. forward increases appreciably the amplitude and the required 
decay time of the initial spin entry oscillation. This effect is 
so great that for the 6 percent c case the steady-state values have 
yet to be attained at the end of 40 seconds (126.4 sec full-scale); 
whereas, for the aft c.g. location the steady state flat spin is 
achieved rapidly and smoothly. 

It has been demonstrated experimentally and analytically 
many times that inertia changes which significantly alter the gyro- 
scopic moments will appreciably modify the spin and recovery charac- 
teristics to be realized with a given aerodynamic configuration. 

The error in measuring the inertia about one axis that would be 
required to produce this effect, however, would not go undetected. 
The errors to be incurred in measuring or estimating the inertias 
would in fact be even considerably less than the differences in- 
curred in attempting .to dynamically simulate a set of full-scale 
values on a dimensionally scaled model. Since the weight and 
inertia must be reduced to the third and fifth power of the scale, 
respectively it is indeed difficult not to exceed the required 
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(desired) values. For example, the I I and I values for the sub- 

Ay z 

Ject model would have to.be reduced by 10.6, 16. 7 and 18.0 percent, 
respectively, to simulate the full-scale airplane. Even though 
these representative scaling errors greatly exceed experimental 
measuring errors, their influence on the spinning motion is pre- 
sented in figure 6. There is a detectable difference between the 
motions. These differences, however, would be considered acceptable 
if obtained during an analytical-experimental correlation study. 

It is interesting to note that since slightly higher rotation rates 
are realized for the lower (desired) inertia, the model results 
are a shade unconservative. 

Significance of Errors in Assumed Aerodynamic Characteristics 

The influence of the static and the dynamic aerodynamic 
characteristics on the spinning motion of an aircraft presented 
in reference 5 were verified during this study. However, 
besides investigating the effect of an individual aerodynamic 
quantity it was also the intent of this study to determine 
the significance on the spinning motion employing two conven- 
tional aerodynamic models which differed in many respects . Both 
models were assumed to represent the same aircraft. Figure 7 
presents time histories which illustrate results deemed to be 
significant and of interest. 

It was found that an equilibrium spin could not be com- 
puted with model 2 since the solutions "blow-up". That is, as 
shown in figure 7a, an ever increasing yaw rate results in a diver- 
gence in (C and y3 . The peak cc value increasingly exceeds 90 degrees 
for which there was no aerodynamic data. Consequently, the time 
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histories computed become invalid and meaningless. As shown in 
figure 7b, the same type of problem is experienced when the rotary 
derivatives based on the low-frequency forced oscillation tests 
are substituted for the model 2 values (high-frequency values). 

A steady flat spin was obtained when model 1 was used 
(see figure 7c). However, when the model 1 rotary derivatives 
were replaeed with the- model 2 values (figure 7d) it was again 
impossible to compute a spin. The same was true when the low fre- 
quency rotary derivatives (figure 7e) were employed in model 1. 

In Appendix C the rotary derivatives used in model 1 are 
compared with the model 2 values in figures C-19 thru C-22. 

Although the derivatives for model 1 and 2 are referred to differ- 
ent moment centers they still are representative of the values to 
be realized when they are all transferred to 14 percent c during 
the spin computation. It can be seen that the values of C n and 
c X r (figure C-21c and C-22b, respectively ) are appreciably different 
for models 1 and 2 in that the model 2 values are unstable at some 
angles of attack.. The importance of each derivative on the motion 
was therefore determined. 

When the model 2 C nr value was employed in model 1, the 
previously computed steady flat spin (figure 7c) was replaced with 
an extremely high frequency oscillatory motion in which the value 
of r is ever increasing with time (figure 7f). This computed motion 
obviously does not represent the "real world". The steady flat 
spin obtained with model 1 was also radically changed when the model 
2 value of C/ r was used in model 1. As shown in figure 7g, a low 
frequency divergent oscillation in a and y3 is realized which event- 
'ually kicks out of the spinning motion. The inability to compute 
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a spin with model 2 was due, therefore, to both C„ and C* 

3 n r x r 

assuming unstable values, over small ranges of angle of attack. 

As shown in figure C-21a, the C m value for model 2 assumed 

Q 

a near zero value above (£,= 75°. When this value was used in 
model 1, the motion presented in figure 7h was obtained. It can be 
seen that this derivative had a significant effect in that the 
steady flat spin became oscillatory and assumed a lower spin rate. 

The C n {Ct i j3) data of model 1 and 2 (see figures C-8 and 
C-l8, respectively) are radically different above an angle of attack 
of 40 degrees. (Both sets of data are presented for a moment re- 
ference center of 14 percent c.) For model 1, the highest level of 
instability (pro-spin yawing moment) attained in the vicinity of 
40 degrees angle of attack is maintained up to 90 degrees. Whereas, 
the pro-spin yawing moment of model 2 rapidly falls off above an 
angle of attack of 40 degrees to a near zero (either slightly pro- 
or anti-spin) value. The spin obtained when the model 2 C n data is 
used in model 1 is shown in figure 7i. The influence these differ- 
ent C n characteristics have on the spin can be seen by comparing 

figure 7i with figure 7c. By effectively removing the large pro- 

o 

spin yawing moment above a = 40 the initial spin entry oscillation 
is more rapidly damped but the developed steady flat spin rate is 
only slightly reduced from 8.7 to 8.2 radians/sec. The reason the 
steady flat spin characteristics were not significantly affected by 
this large aerodynamic change is because the aerodynamic yawing moment 
approaches a value clo8< to zero when spin equilibrium is achieved. 
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It was concluded from this study that the same type of 
spin could be obtained with the conventional aerodynamic models 
1 and 2 if the stable Cj^, and C m ^ derivatives of model 1 
were used in model 2. 

Significance of Selecting a Different Type of Analytical Model 

As stated previously, the conventional technique does not 
include the aerodynamic forces and moments acting on a spinning 
aircraft due to steady rotational flow nor does it limit the con- 
tribution of the rotary derivatives to the oscillatory component 
of the total angular rates. Whereas, the analytical technique 
employing rotation— balance data does both. It was desired, there- 
fore, to enter a spin using the same controls employing both tech- 
niques. Because of the limited rotation-balance data available, 
this desire was achieved using the procedure previously discussed 
and by limiting the correlation to the developed spin and recovery 
characteristics. 

Since the available rotation-balance data was limited to a 
positive sideslip angle, the initial attempt was to calculate left 
spins which would assume small positive values of sideslip. Equili- 
brium spins, however, could not be obtained with aerodynamic model 
3 °r When only the /i = 0° data was employed with model 4, 

® ^*l^t spin was attained, but this spin also could not be 
maintained. An analysis of the rotation-balance data 
indicated anomolies in the counter-clockwise data (left spin) 
for both = o and 5 degrees which could not be explained. 

It was decided, therefore, to compute right spins using only 
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y3= 0 rotation-balance data. 

Figures 8a and b present the time histories for the motions 
attained using the conventional technique (model 1) and the technique 
based on rotation-balance data (model 3), respectively. (Model 3 
consists of the rotation-balance data and the rotary derivatives 
of model 1.) In both cases, a steady flat spin is obtained. How- 
ever, a more severe flat spin is obtained with the rotation-balance 
set-up. The angle of attack being increased from 76.7 to 85.4 deg- 
rees and the spin rate from 5.4 to 9.7 radians/sec. 

Figure 8c presents the time history computed with the rota- 
tion-balance data and the rotary derivatives from model 2; this data 
package being referred to as model 4. With model 4, the flat span 
obtained became oscillatory and spun at a lower angular rate. How- 
ever, as discussed in the preva’ous section, it was impossible to 
compute a spin with this set of rotary data when the conventional 
analytical technique was used. It should also be noted that in 
this instance, substituting the model 3 C nr or Or value in model 4 
had absolutely no effect on the motion. The effect of using the 
Ciriq v al' ie of model 3 in model 4 can be seen by comparing figure 8d 
with 8c. With this one change, the spins obtained with models 3 and 
4 are almost the same. The angle of attack and spin rate are 84.8 deg 
and 9.14 radians/sec, respectively. The full-scale values would 
correspondingly be 84.8 deg and .46 revolutions/sec which compares 
very favorably with the 87 deg and .47 revolutions/sec values ob- 
tained with the i - scale spin model. These results were obtained 

36 

for almost the same i s , 5 a , 6 R control configuration. The only 

difference being that the spin model was tested with 6 = + 7° and 

a 
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the rotation-balance data was obtained for & =4-5°, When both 

the C /p and C np values were replaced in model 4 with those from 
model 3 , &n insignificant effect was noted (compare figure 8d with 
&e ) . 

These results 

1) strongly support the belief that the conventional analy- 
tical technique employed to date is a pseudo-technique 
which cannot be used to predict full-scale spins and, 

2) indicate that past conclusions (reference 5) analyti- 
cally arrived at relative to the significance of rotary 
derivatives on the spinning motion are in many instances 
wrong. 


A direct comparison between the recovery characteristics 
using the conventional and rotation-balance data technique is not 
possible since a more severe (higher a and .a) flat spin was obtained 
with the rotation-balance data. However, from figure 9 it does not 
appear that any significant difference is realized relative to the 
effectiveness of the recovery controls. In both cases, the recovery 
controls were introduced at t = 20 seconds and it took six seconds 
to stop the spin using the conventional technique. An additional 
two seconds was required to terminate the more severe rotation- 
balance data spin. 


A Simple Aerodynamic Model Based on Rotation-Balance Data 


A preliminary investigation indicated that it is possible 
to extract static aerodynamic derivatives (i.e., C m , C, , C 


'n 6 , 


» ^m^ 


s 


l a v 3 

etc.) from rotation-balance data. If this is indeed 
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2h 
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the case for all aircraft configurations, then 

1) rotation-balance tests eliminate the need for corres- 
ponding static wind tunnel tests and 

2) a simple aerodynamic model can be constructed which is 
more valid than the conventional aerodynamic model since 
it includes the aerodynamics associated with rotational 
flow and, as discussed previously, permits one to apply 
the rotary derivatives in a more proper manner. 

With a simple model one can employ all of the unstalled and stalled 
static ( — o) aerodynamic data that may be available and re- 

moves the issue of how to smoothly transition from one type of aero- 
dynamic model to another during the incipient spin calculation. 

Also, the model can easily be incorporated into existing large- 
angle six-degree-of-f reedom computer programs. 

The following simple model would adequately represent the 
vehicle investigated herein. 


'm* 


as a function of cc 
as a function of cc,/3 


C X , C n , Cy as a function of cc,/3 

c n> c y as a function of , <x, i s 

C m< , C c , C N as a function of Cc 
1 s is is 

c / t c . c n< - , > c y. „ as a function of CC . -CL b 
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CONCLUSIONS 

The following major conclusions are based on the results 
of the analytical studies reported herein. With the exception of 
the first, all of the conclusions pertain to factors which can 
influence the outcome of an experimental-theoretical correlation 
study. Although these studies were made for one configuration, 
the conclusions should be applicable to other aircraft. 

1. Experimental spin investigations conducted in a con- 
stant density environment duplicate the Froude number only at the 
initial full-scale spin altitude since the full-scale airplane at 
high altitude experiences large atmospheric density changes during 
the spin. Therefore, the full-scale rate of sink and spin rate 
predicted from a constant density environment test will be greater 
(conservative) than the actual airplane values. 

2. The spinning motion can be acutely sensitive to a small 
change in the application of the spin entry controls. Since the 
critical control phasing required to trigger a radical change in 
the spin behavior of the aircraft cannot be predicted at this time, 
flight recorded control time histories must be faithfully duplicated 
during an experimental-analytical correlation study. 

3. The spinning motion, both the incipient and developed 
phases, are sensitive to the location of the center of gravity but 
are not significantly, effected by errors incurred in measuring 
airplane inertias. 

4. The presence of a bank angle at the time spin entry con 
trols are applied does not influence the devleoped spin. The inci- 
pient spin motion, however, is affected by bank angles 40 degrees 
and greater in that entry into the developed spin is achieved more 





smoothly if the aircraft is initially banked in the direction of 
the spin. 

5. The damping in pitch derivative, c rriq> has a signi- 
ficant influence on all phases of the spinning motion. The magni- 
tude of the static yawing moment in the 50 to 90 degree angle- of- 
attack range appreciably affects the incipient spin but only 
slightly affects the developed spin. 

6. A spin may not be computed using the conventional analy- 
tical technique when the rotary derivative C m , C„ or C a is un- 

m q n r " r 

stable over some portion of the angle -of-attack region. This is not 
the case, when the contribution of the rotary derivatives are limited 
to the oscillatory component of the total angular rates which is the 
proper procedure. Some of the conclusions arrived at in the past 
usi-Hc the conventional analytical technique, therefore, are incorrect 
relative to the significance of rotary derivatives on the spinning 
motion. 

7. The spin computed with rotation-balance data duplicated 
the developed spin obtained in the spin tunnel. The spin computed 
in the conventional manner with static aerodynamic data was of the 
same type (flat) but was considerably less severe. It is indicated, 
therefore, that the aerodynamic moments generated in the spin due to 
steady rotational flow, as measured by a rotation-balance, are in- 
deed significant. The analytical technique used to date is, there- 
fore, a pseudo-technique which cannot be used to predict full-scale 
spin motions. 
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8. It appears that static aerodynamic derivatives can 
be extracted from rotation-balance data. It is possible, there- 
fore, to simply add the aerodynamics associated with steady rota- 
tion flow to the conventional aerodynamic model. In this man- 
ner, not only is the aerodynamic model more valid, but. i.1, also 
permits one to modify the mathematical model to more properly 
apply the rotary derivatives. 


C 
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Figure 2.- Influence of rate of change in atmospheric density with altitude during full- 
scale spin. 
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Figure 4,- Sensitivity of spinning motion to initial Euler angles. 
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Figure 5.- Sensitivity of spinning motion to center-of-gravity location. 
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a) 6 percent c 
Figure 5.- Continued. 
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Figure 5.- Continued. 
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Figure 6.- Concluded, 
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a) Model # 2 


Figure 7.- Influence of some aerodynamic characteristics on spinning motion. 
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a) Model # 2 
Figure 7.- Continued, 
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b) Model # 2 with low frequency dynamic derivatives. 


Figure 7.- Continued. 
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i Figure 7,- Continued. 
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•) Model # 1 with low frequency dynamic derivatives. 


Figure 7,- Continued. 


G\ 


OMEGA 



BETfl-OEG ALPHR-OEG FLIGHT PATH ANGLE 



f) Model # 1 with high frequency (Model # 2) C n 
Pigure 7 , - Continued, 
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g) Model # 1 with high frequency (Model # 2) Cj 


Figure 7.- Continued. 
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g) Model # 1 with high frequency (Model # 2) 
Figure 7,- Continued. 


SI 

ro 



h) Model # 1 with high frequency (Model # 2) C m . 

q 


Figure 7.- Continued. 
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Figure 7.- Continued. 
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i) Model # 1 with Model |2C n 
Figure 7.- Concluded. 
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Figure 8.- Continued, 


CD 





•si 

VO 















RLPHfl-DEG FLIGHT PRTH RNGLE 



Figur* 8,- Continued, 
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Figure 8.- Concluded 
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APPENDIX A 

EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 


The dynamical equations required to specify the transla- 
tional and rotational motions of a rigid body moving through space 
are described in this appendix. The familiar six degree of freedom 
differential equations representing linear and angular accelerations 
of a moving body axis system having its origin at the aircraft center 
of mass are presented below. 
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In addition, the following formulas were used : 
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Turns in spin 
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dt 
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©e “ q cos 0 e - r sin 0 e 

0 e = p + r tan e e cos 0 e + q tan © e sin 0 € 
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* * 3r + 9 0 


r + r 
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These total angular velocities (p, q, r) consist of steady 
rotation (p r , q , r r ) components upon which oscillatory (p q and 
r Q ) components are superimposed. These components are def£ned°as 
follows : 


P r * -Ye sin © e 

<Jr = Y e cos © e sin 0 e 
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For the conventional model, the following total 
derivatives were used 
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For the rotation balance model, the following total 
derivatives were employed, where C , C , C , etc. are rotation 
balance data -*• y z 
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SYMBOLS 
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Measurements and calculations were made in the U.S. Custo- 
mary Units. Factors for converting these units to the International 
Systems of Units may be found in reference 6 . 

The Body Axes System is used with the origin of the axes 
system located at the aircraft center of gravity. The X axis is paral- 
el to the aircraft Fuselage Reference Line (FRL) and is positive forward, 
the Y axis is positive towards the right wing tip and the Z axis is 
positive downward. The following illustration shows this axis system 
and the positive direction of the angles, forces and moments associated 
with it. 



In the following definitions a dot (*) over a quantity de- 
notes one differentiation with respect to time and a zero subscript ( ) 
denotes the initial value of the quantity. ° 
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c 

1 x» I y» I e 

moments of inertia about the X,Y and 2 
body axes, respectively 

slug-ft 2 


1X2 

product of inertia, positive when the 
principle X axis is inclined below the 
X body axis at the aircraft nose 

slug— ft 2 


m 

aircraft mass (=W/g) 

slug's 


w 

aircraft weight (=mg) 

lbs 


s 

wing area 

OJ 

■p 


b 

wing span 

ft 


c 

wing mean aerodynamic chord 

i 

ft 


w/s 

wing loading 

lb/ft 2 

c 

(Iy - I z )/mb 

inertia rolling moment parameter 



( x z ~ I x )/mb 2 

inertia pitching moment parameter 



( X x ” Iy)/mb 2 

inertia yawing moment parameter 

- 


t 

time 

sec 


1° 

q 

atmospheric density 

2 

pY 

dynamic pressure (~ 

slug/ft 3 

lb/ft 2 


H 

vertical height 

ft 

c 

<£ 

angle of attack, measured between the X 
body axis. and the projection in the X-Z 
plane of the relative wind vector, positive 
when the X axis is above the projected 
relative wind vector 

deg 
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fi 

angle of sideslip, measured between the 
relative wind vector and its projection 
in the X-Z plane, positive when the re- 
lative wind vector is to the right of the 
X-Z plane 

deg 

®e 

inertial pitch attitude angle, measured in 
a vertical plane between the X body axis and 
the horizontal plane, positive when the X 
axis is above the horizontal plane 

deg 


inertial roll attitude angle, measured in the 
Y-Z plane between the Y body axis and the 
horizontal plane, positive when the Y axis 
is moving clockwise as viewed from behind 
the aircraft 

deg 

♦'e 

inertial yaw attitude angle, measured between 
the initial flight direction and the projec- 
tion in the horizontal plane of the X body 
axis, positive when the X axis is moving 
clockwise as viewed from above 

deg 

0 

angle between Y body axis and horizontal mea- 
sured in vertical plane, positive for erect 
spins when right wing is downward and for 
inverted spins when left wing is downward 

deg 

r 

flight path angle, measured in a vertical plane 
between the horizontal plane and the resultant 
velocity vector, positive when the resultant 
velocity vector points above the horizontal 
plane 

deg 


differential horizontal tail deflection (half 
angle), positive to produce left rolling moment 

deg 


control deflection-stabilizer, positive direc- 
tion trailing edge down 

deg 

*R 

control deflection-rudder, positive direction 
trailing edge left 

deg 


C 
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v a 

v 

U, V, w 


u 


> 



V, q, r 


total linear velocity vector 
free-strearn velocity 

components of the total linear velocity 
vector (V R ) along the X, Y and Z body 
axes, respectively 


components of the relative linear 
acceleration vector along the X,Y 
and Z body axes, respectively 


components of the total angular velocity 
vector along the X, Y and Z body axes, 
respectively 


frequency of for ed-oscillation tests 




total angular velocity vector 


g 


acceleration due to gravity 


Fx ae ro * Fy ae ro * 

F z 

aero 


components of the aerodynamic forces along 
the X, Y and Z body axes, respectively, 
excluding the contributions due to the 
rate derivatives Cn- and C v , 

* y /S 


L , M , 
aero aero 

N aero 


components of the aerodynamic moments 
about the X, Y and Z body axes, 
respectively 


T 



-F x 

aero 

q S 


total engine thrust force 

total chordwise force coefficient, 
positive along the negative X body axis 


ft/sec 
ft/sec 
ft /sec 

ft/sec 2 

radians/sec 

radians/sec 

radians/sec 

ft /sec 2 

lbs 

ft- lbs 
lbs 
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c = . 

aero 

y 

q S 

/ 

-f 2 

c lC 

aero 

q S 

C> 

L 

aero 


QSb 

<- 

M 

aero 

qSc 


total sideforce coefficient, excluding the contribution 
due to C y ^. , positive along the positive Y body axis 


total normal force coefficient, excluding the contribu- 
tion due to positive along the negative Z body axis 


total rolling moment coefficient, positive direction of 
moment drives the right wing tip down 


total pitching moment coefficient about the aerodynamic 
center, positive direction of moment drives the nose up 


/ N 

C n = — - a - ero total yawing moment coefficient, positive direction of 

qSb moment drives the nose right 


C c (ot) 

C 

C y (y3,a;) 


variation in chordwise force coefficient due to ancle 
of attack e 


variation in side force coefficient due to sideslip 
angle and angle of attack 


C N (cc) 


variation in normal force coefficient due to angle of 
attack 


C Z (/$>&) 


variation in rolling moment coefficient due to sideslip 
angle and angle of attack 


Cm (*) 


variation in pitching moment coefficient due to ancle 
of attack 


c n (y3,o:) 


variation in yawing moment coefficient due to sideslip 
angle and angle of attack 



C v = 




C. = 


-58 



per deg . 


per deg 
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3J3 




per deg 


per radian 


per radian 



per radian 


per radian 


per radian 


per radian 


per radian 


per radian 


C 
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° y *a 

s 

Jfz_ 

»»a 

per deg 

C ^a 

= 

i°jC 

J 6 a 

per deg 

n 


ac y 

per deg 

n «a 


^a 

c 


ac 

y 

per deg 

y *R 



c. 


ac. 

per deg 



d6 R 

C„ 


* C „ 

per deg 

6 R 


j 6r 

Cr 


ac c 

per deg 

C i S 

■ - 



C N4 

X S 


a c 

per deg 


b ^ s 

0, , 

s 

a 0 . 

per deg 

m 1 . 
• • 


ai (5 

c 

7 

- 

- C N 

per deg 

,: X 
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Figure C-4. - Effect of angle of attack on lateral-control effectiveness. Lamda 
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Figure C-5. - Continued. 





Figure C-5. - Concluded 

























































Figure C-9. - Effect of angle of attack on damping derivatives 









Figure C-9. - Continued 
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Figure C-10. - Concluded 
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Figure C-ll. - Aerodynamic derivatives assigned zero values. 
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Figure C-15 
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C-17. - Effect of angle of attack and longitudinal- 
deflection on rolling-moment coefficient. Lamda = 22°. 
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Figure C-21 . - Effect of angle of attack and forced-oscillation frequency on damping derivatives. 
Lamda = 22°. 
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a) 0 = 0 , e = 55, 6o, 6 5 

Figure C-24. - Effect of rotation rate and pitch and roll attitude angle 

on the longitudinal coefficients for i s = -30°. Lamda =22 
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Figure C-25. - 


a) 0m o“, e = 55, 60, 65° 

n5 f ?h» ?£ rotation rate and pitch and roll attitude an^le- 
t?ols? 2da"= 2|° coefflclents f°r right pro-spin con- ‘ 










■ H 

■■■ 


o 70 


□ 75 


O 80 




■MM 


BSIlSIsI*! 









0, deg 


O 55 


□ 60 


O 65 







»S 


issSSSSS» 


Figure C-26. - 


a) 0 = 0 = 0°, e = 55, 60, 65° 

Effect of rotation rate and pitch and roll attitude 
angles on the longitudinal coefficients for left pro 
spin controls. Lamda = 22° * y 
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Figure C-28. - Concluded. 

/ / 


/ 






2 


3 


90 


o 


tinued 










EiSISISs 


□ 90 


I—— m aaM 

TiirPfTi— H—i^a 




mmmm 


SHil 





.2 


3 


and 90° 


tinued 





2 


3 


90 ° 
uded . 




